SCIENTIFIC 

REPpRTS 



»H5 




SUBJECT AREAS: 
CHEMICAL BIOLOGY 
DRUG DISCOVERY 
PHARMACOLOGY 
SENSORS 

Received 
4 May 2011 

Accepted 
23 June 2011 

Published 
7 July 2011 



Correspondence and 
requests for materials 
should be addressed to 
Y.F. (fangy2@corning. 



Label-free integrative pharmacology 
on-target of drugs at the ^-adrenergic 
receptor 

Ann M. Ferrie, Haiyan Sun & Ye Fang 

Biochemical Technologies, Science and Technology Division, Corning Inc., Corning, NY 1 483 1 , USA. 

We describe a label-free integrative pharmacology on-target (iPOT) method to assess the pharmacology of 
drugs at the p 2 -adrenergic receptor. This method combines dynamic mass redistribution (DMR) assays 
using an array of probe molecule-hijacked cells with similarity analysis. The whole cell DMR assays track cell 
system-based, ligand- directed, and kinetics-dependent biased activities of the drugs, and translates their 
on-target pharmacology into numerical descriptors which are subject to similarity analysis. We demonstrate 
that the approach establishes an effective link between the label-free pharmacology and in vivo therapeutic 
indications of drugs. 

Embraced by the "one gene, one drug, one disease" philosophy, drug development campaigns heavily invest 
on high throughput screening (HTS) -compatible molecular characterization assays to search maximally 
selective drugs ("magic bullets") for target-based therapies 1 . However, the best in class drugs are identified 
more by serendipity than by rational design 2 , and prioritization of lead candidates is as much an art as a process 3 . 
This is mostly due to the poor correlation between in vitro results and in vivo indication(s). To prioritize lead 
candidates and predict their therapeutic potentials, biological fingerprints based on binding profiles 4 , gene 
expression profiles 5 ' 6 , cellular phenotypic effects 7 , side-effects 8 , and chemical structures 910 can be produced so 
that these molecules can be directly compared and sorted. These fingerprint-based approaches are promising to 
relate molecular, cellular and in vivo features of drug molecules to targets and diseases. However, these 
approaches generally have poor resolution in the on-target pharmacology of a drug, the functional consequences 
of the drug binding to a specific target. The binding profile of a drug does not directly translate to its functional 
activity. The cellular phenotypic responses can differentiate antagonists from agonists, but suffer significantly in 
target specificity, in part due to the presence of compensatory pathways and network interactions, and in part due 
to the accuracy of the reference signature associated with a specific target or a cellular process 5 . Furthermore, 
diverse molecular assays have revealed pluridimensional efficacy {i.e., an assay readout dependent efficacy) of a 
single drug acting through a specific target 1112 . Since these assays measure single end points downstream the 
target activation, it is obviously challenging to directly compare and prioritize molecules based on multiple 
individual assay results, thus making it difficult to comprehend the therapeutic potential of these molecules. 
These problems may be overcome if it were possible to effectively cluster drugs based on their in vitro multifaceted 
on-target pharmacology. 

G protein -coupled receptors (GPCRs) represent the largest and most successful class of druggable targets in 
human genome. This is illustrated by (3-adrenergic receptors ((3-ARs). The (3-ARs belong to class A GPCRs and 
consists of Px-, p 2 - and P 3 -AR subtypes. Drug development against P-ARs has been fruitful in the past several 
decades (Supplementary Table SI). P-blockers have been used for the treatment and management of cardio- 
vascular conditions 13 , migraine 14 , and ophthalmic disorders 15 , presumably due to blockage of the activation of 
Px -receptor 16 . p 2 -agonists have been long used in and still are the most effective bronchodilators for the treatment 
of asthma 17 . Several P-agonists are also used for the treatment of cardiac decompensation, anaphylaxis, sepsis, and 
premature labor 18 . However, neither all beta-blockers behave equally for treating various heart diseases, nor are all 
P 2 -agonists effective in the management of asthma. Together with diverse in vivo indications, the distinct clinical 
profiles of p-drugs have challenged our view of how drugs should be classified, which, in turn, should guide us 
how to screen and test drugs in vitro 19 . 

Here we present a strategy to determine on-target pharmacology of drugs acting on GPCRs using label-free 
integrative pharmacology on-target (iPOT). Central to iPOT is label-free whole cell dynamic mass redistribution 
(DMR) assay 20-23 . DMR assays use a label-free optical biosensor, as shown in Fig. la, to non-invasively record a 
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Broadband light Reflected light 

Figure 1 | The principle of the iPOT. The iPOT combines DMR assays with similarity analysis to assess on-target pharmacology of the (3-AR drugs, 
(a) Resonant waveguide grating biosensor for receptor signaling. The biosensor uses leaky mode nano-grating waveguide structure to generate an 
evanescent wave to sensor whole cell responses. Cells are directly cultured onto and become adherent to the biosensor surface via adhesion complexes. A 
schematic of (3 2 -AR signaling pathway is also included. The intervention with various probe molecules of distinct signaling proteins in the (3 2 -AR 
pathway can be used as the basis to determine the on-target pharmacology of various |3-drugs. (b) Assay protocols that use various probe molecules to 
precondition a cell expressing the (3 2 -AR for reporting DMR pharmacology of the (3-drugs. (c) A numerical descriptor of a (3-drug-induced DMR. The 
responses at four distinct time points were extracted and color coded. Red: positive value; green: negative value; black: close to zero. The DMR was 
obtained by stimulating native A431 with salmeterol at 10 uM. 



drug induced dynamic redistribution of cellular matter within 
— 150 nm of the sensor surface 20 . The DMR is a very discerning assay 
that gives highly texture data at the whole cell level 19,22 . A drug- 
induced DMR is often sensitive to cellular context, which can be 
controlled using genetic manipulations or interventions with small 
probe molecules. DMR measurements in these perturbed cells can 
manifest the sensitivity of a drug-induced DMR to the target(s) /path- 
way^) with which the small probe molecules intervene. Combining 
DMR assays in a wide array of probe molecule perturbed cells with 
similarity analysis can separate drugs based on their distinct on- 
target pharmacology. We demonstrate that this strategy provides 
an effective link between in vitro profiles and in vivo indications of 
P- drugs. 

Results 

The principle of iPOT. The iPOT begins with the use of diverse 
probe molecules to hijack a cell or cell system, followed by 
profiling drugs with DMR assays. The DMR profiles obtained for 
each drug are then translated to a multi-dimensional coordinate such 
that all drugs tested can be compared using similarity analysis. The 
probes are chosen to recapitulate signaling pathways downstream 
the target as well as pharmacological activities of drugs, such that 
the pathway biased activity, if any, of drugs can be systematically 
surveyed. The probes can be toxins for G proteins, inhibitors for 
kinases and activators for enzymes within the receptor signaling 
cascades, or P-drugs themselves (Fig. la). The DMR arising from 
P-drugs in the probe pretreated cells are used directly for similarity 
analysis, except that when the probes are P-drugs, the DMR arising 
from epinephrine or propranolol are used as readouts (Fig. lb; 
Table 1). The cell or system is chosen based on the known 
signaling capacity of the target, although the cells that are derived 



from in vivo action sites of the drugs can be used. Hijacking of the 
cell with the probes redirects the signaling routes of the target, 
thus manifesting the biased activity of drugs towards the probe- 
intervened pathways. After translating the kinetic DMR profiles 
into multidimensional coordinates, similarity is analyzed to 
categorize drugs into distinct clusters. We assume that drugs 
within a given (sub)cluster share a common mode of action 
possibly linked to an in vivo indication. The iPOT effectively 
integrates the system-based, ligand-directed, and kinetics- 
dependent biased activities of the drugs. 

iPOT profiles established a link between in vitro results and in 
vivo indications of p-drugs. We implemented the iPOT to 
characterize all P-drugs approved by the Food and Drug 
Administration (FDA) (Supplementary Table S2). All p 2 -drugs 
were profiled at 10 |iM to achieve maximal signaling capacity and 
to be amenable for HTS (Supplementary Table S3). Human 
epidermoid carcinoma A431 was chosen to be the cell system since 
its endogenous P 2 AR has been served as a model to elucidate the 
P 2 AR signaling in the past 24 . Before measuring the DMR induced by 
p-drugs, A431 was preconditioned under eight different conditions: 
(1) the assay vehicle for 1 hr as a control; (2) 1 |iM propranolol (PRO) 
for 1 hr; (3) 5 nM epinephrine (EPI) for 1 hr; (4) 5 nM EPI for 5 hrs; 
(5) 400 ng/ml Chorea toxin (CTx) for overnight; (6) 100 ng/ml 
Pertussis toxin (PTx) for overnight; (7) 10 |iM 4,5,6,7-tetrabromo- 
benzotriazole (TBB) for 1 hr; and (8) 10 |iM U0126 for 1 hr. We also 
measured the EPI DMR in cells after pretreated with 10 |iM P-drugs 
for 1 hr or 5 hrs, the p-drug DMR in the presence of 10 |iM forskolin 
(FSK), and finally the PRO DMR in cells after pretreated with 10 |aM 
P-drugs for 5 hrs (Fig. lb). Total twelve DMR profiles were collected 
for each drug (Table 1). 
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Table 


1 Assay protocols and DMR signals used for 


similarity and correlation anc 


lysis. 


Assay 


number Probe, pretreatment duration 


DMR readout 


Note 


1 
i 


0 1 °L DKA^d 1 hr 


Ligana, 1 U ji/V\ 


rositive control, agonism 


o 
z 


1 0 |iM Licjcind, 1 hr 


Epinephrine, 5 nM 


Specificity, relative potency and efficacy, modes of action 


Q 

o 


i |_i/v\ propruiKjicji, 1 111 


1 inriirifA 1 M i i KA 

Liyuiiu, i \j \iiv\ 


OpcLIIILITy, IcIUTIVc pOTcllLy UlIU cMILULy, iTlvJUco (Jl ULTIvJIl 


4 


U ll/Vi E Ul NcUll 1 1 lie, 1 III 


1 innnrl 1 i i 1 1 KA 
LIUUIIU, 1 \J \Xly\ 


IXcCcpiUI Ucbcl lol 1 IZ.U 1 1 (Jl 1 UlIU UcUtl 1 VU 1 1 Ul 1 


o 


1 0 |iM ligand, 5hr 


Epinephrine, 5nM 


Kinetics dependent receptor blockage by antagonists and 








receptor desensitization by agonists 


A 
U 


j ri/v\ tpiricprir iiic, vjiir 


1 1 n m r\ i n 1 1 ka 
Liyuiiu, i \j \iiv\ 


IxlllcTILo UcpcriUcllT UcULTI VUTIOll U\ Tile IcLcpTUI Unci 








piUIUIIUcU ULIIVUIIUII Willi epi 1 lepi II 1 1 Ic 


7 


1 U |i/V\ ligana, onr 


d l l in 
rropranolol, 1 U |i/V\ 


Relative potency, kinetics dependent agonism 


8 


0.1% DMSO, lhr 


Forskolin 1 0 jiM 


cAMP-PKA pathway 






Ligand, 1 0 jiM 


9 


400 ng/ ml cholera toxin, 20hr 


Ligand, 1 0 |iM 


G as -dependent and independent activity 


10 


1 00 ng/ ml pertussis toxin, 20hr 


Ligand, 1 0 |iM 


G ai sensitivity 


11 


lOuMTBB, lhr 


Ligand, 1 0 jiM 


CK2 sensitivity 


12 


IOjiM U0126, lhr 


Ligand, 1 0 |iM 


MEK1/2 sensitivity 



Each DMR profile is a kinetics response with more than 200 time 
points. To improve the efficiency of similarity analysis, we reduced 
the DMR dimensions to four distinct time points (3, 5, 9, and 45 min 
post stimulation) (Fig.lc). The dimension reduction is based on 
clustering of time domains of the DMR of all p- drugs in the 
DMSO treated cells (Supplementary Fig.Sl). Results showed that 
the p-drug DMR generally propagate with four distinct time periods: 
immediate (2-4 min), early (5-6 min, and 7-9 min), and late res- 
ponses (10-60min post stimulation). Similar pattern was obtained 
for all DMR signals under all twelve conditions. Thus, we selected 
one time point from each period to represent each DMR, and we 
found that these time points adequately represent the key features of 
the drug responses (Supplementary Fig.S2 and S3). Thus, the twelve 
DMR profiles of each drug can be rewritten into a 48 dimensional 
coordinate. 

Similarity analysis using unsupervised Ward hierarchical cluster- 
ing algorithm and Euclidean distance metrics 25 categorized drugs 
into distinct clusters (Fig.2). Detailed analysis revealed several inter- 
esting aspects. First, the DMR assays are high quality in general, since 
two isomers and their mixture of propranolol led to subtle difference 
in DMR profiles but were still bucketed together; and the same was 
true for S-isomer and its mixture of pindolol. Further, the p 2 -inactive 
R-( + ) -atenolol was clearly different from its active isomer, S(-) -ate- 
nolol. R-(+)-atenolol was clustered with practolol and the negative 
control dimethyl sulfoxide (DMSO). Practolol is a Pi-selective but 
much less potent p 2 -antagonist. Second, two superclusters were 
formed; one for partial agonists and antagonists, and another for 
strong partial agonists, full agonists and the adenylyl cyclase activator 
forskolin (FSK). Each supercluster is further divided into many sub- 
clusters. Third, drugs in each subcluster share almost identical thera- 
peutic indications as approved by the FDA and documented in 
DrugBank 26 , and the correlation between DMR profiles and in vivo 
indications was correctly established for more than 90% of the drugs 
tested. Notable disassociations between the DMR profiles and ori- 
ginal indications were: (1) carvedilol was clustered with PRO. A small 
but positive clinical trial using carvedilol for the prophylactic treat- 
ment of migraine has been reported in literature 27 . (2) Terbutaline 
was close to ritodrine, and to less degree, phenylephrine. Although it 
is controversial and has not been approved by the FDA for tocolytic 
use, terbutaline is the most widely prescribed beta-mimetic to arrest 
or prevent premature labor in the United States 18 . (3) Salmeterol was 
distinct from all other anti-asthma agonists. (4) Timolol was distinct 
from PRO although timolol is also used for treating migraine. (5) The 
active metabolite (methylnorepinephrine) of methyldopa was close 
to norepinephrine. 

We further compared the results obtained with DMR profiles with 
those obtained using in vivo indications or chemical features. 



Clustering of P-drugs based on their in vivo indications, as documen- 
ted in DrugBank, led to a relatively low resolution heat map (Fig.3a). 
For this analysis, a binary numeral system, 1 for an approved indica- 
tion and 0 for none, is used to describe each drug based on their 
clinical indications. We also analyzed the chemical similarity of all 
P-drugs tested using ChemMine 28 . ChemMine is an online service for 
analyzing and clustering small molecules. Chemical similarity based 
on Tanimoto coefficient led to bucketing of compounds that shows 
relatively poor correlation with their in vivo indications (Fig.3b). 
These comparative analyses suggest that the iPOT gave rise to better 
correlation between in vitro results and in vivo indications. 

Correlation analysis identified the origin of similarity among 
p-drugs. We used correlation analysis between distinct conditions 
to further determine the origin of similarity/difference among 
P-drugs. We used seven assays to separate P-drugs based on 
agonism, specificity, potency, and kinetics -dependent mechanisms 
of activation and inactivation. First, we recorded the DMR arising 
from all drugs in the native cells. Results showed that there were three 
classes of drugs: full agonists that triggered a DMR similar to the EPI 
DMR; partial agonists that led to a DMR smaller than the EPI DMR; 
and apparent inactive drugs that led to a net- zero DMR 
(Supplementary Fig.S3). Second, we examined the impacts of P- 
drugs on the succeeding EPI DMR, wherein the cells were 
pretreated with P-drugs for lhr or 5hr. Results suggest that S-(-)- 
atenolol, (±)atenolol and practolol exhibited weak antagonism, 
while dobutamine and dopamine led to weak agonism, but R-( + )- 
atenolol was inactive (Fig.4a). However, the EPI DMR after the cells 
were pretreated with P-drugs for 5hr showed that dobutamine and 
dopamine caused a complete desensitization, but bisprolol and 
acebutolol led to a reduced blockage of the succeeding EPI DMR 
(Fig.4a; supplementary Fig.S4). Third, we examined the DMR 
induced by P-drugs after the cells were pre-activated with 5 nM 
EPI for 1 hr or 5hr. Results showed that there was a linear 
correlation between the two DMR for most P-drugs (Fig. 4b). The 
drugs in the cells after the longer EPI pretreatment generally led to a 
more positive DMR. However, the DMR arising from satolol, 
R(+) atenolol, carvedilol, and (±)-, R-( + )- and S-(-) -propranolol 
were largely insensitive to the pretreatment duration with EPI 
(Fig.4b; Supplementary Fig.S5 and S6). Fourth, we examined the 
ability of PRO to block the DMR arising from P-drugs, wherein 
PRO was used to pre- stabilize the receptor at its inactive 
conformation. This further separated p-drugs based on potency 
(Fig.5a). The percentage of inhibition by PRO was mostly inversely 
correlated with the potency (i.e., K { or EC 50 , c amp) of an agonist drug, 
except for a group of partial agonist drugs that include pindolol, 
alprenolol, oxprenolol, labetalol and isoetharine (Supplementary 
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Figure 2 | DMR heat map of clinically available P-drugs. This heat map was obtained using DMR profiling of the drugs under twelve conditions, 
followed by similarity analysis using the Ward hierarchical clustering algorithm and Euclidean distance metrics. The twelve DMR used for analysis were: a 
drug-induced DMR in the cells pretreated with the assay vehicle for 1 hr (Buffer - AR), 1 uM propranolol for 1 hr (PRO - AR), 5 nM EPI for 1 hr 
(EPI - AR), 5 nM EPI for 5 hrs (EPI5h - AR), 400 ng/ml CTx for overnight (CTx - AR), 100 ng/ml PTx for overnight (PTx - AR), 10 uM TBB for 1 hr 
(TBB - AR), and 10 uM U0126 for 1 hr (U-126 - AR). In addition, we also included an EPI-induced DMR in the cells pretreated with 10 uM (3-drugs for 
1 hr (AR - EPI) and for 5 hrs (AR5h - EPI); an PRO-induced DMR in the cells pretreated with 10 uM P-drugs for 5 hrs (AR5h - PRO); and a forskolin 
(FSK)-induced DMR in the presence of (3-drugs (FSK + AR). 



Fig.S7). Finally, we examined the ability of PRO to alter the course of 
a drug-induced DMR after the drug was exposed to the cells for 5 hrs 
(Fig.5b). This further separated agonist drugs, in which PRO 
reversed the DMR arising from distinct agonist drugs to different 
degrees (Supplementary Fig.S8). 

Next, we used three assays to delineate drug pharmacology based 
on G protein -dependent and independent signaling. First, we exam- 
ined the DMR induced by FSK in the absence and presence of a 
P-drug (Fig.6a), wherein FSK was used to activate adenylyl cyclases. 
This delineated drugs into three groups: drugs that suppressed the 
FSK DMR, drugs that potentiated the FSK DMR with an accelerated 
kinetics, and drugs that had little impact on the FSK DMR 
(Supplementary Fig.S9). Notable is salmeterol which significantly 
potentiated the FSK DMR, unlike all other strong partial and full 
agonists. Second, we compared the p-drug DMR without and with 
the pre- decoupling of G as proteins from the receptor by CTx 
(Fig.6b). This identified drugs that still produced noticeable DMR 
(e.g., EPI, alprenolol), or led to an increased DMR (e.g., carvedilol) in 
the CTx-pretreated cells (Supplementary Fig.SlO). Third, we com- 
pared the P-drug DMR without and with the masking of G ai proteins 
by PTx (Fig. 7a). This resulted in a linear correlation between the two 
DMR, but with a slope of 0.95, suggesting that the impairment of G ai 
by PTx suppressed all p-drug-induced DMR. 



Finally, we used three assays to probe ligand- directed biased agon- 
ism. First, we used a MEK1/2 inhibitor U0126 to inhibit the activity 
of MEK1/2. This led to a linear correlation with a slope significantly 
greater than 1 between the two DMR arising from P-drugs without 
and with the U0126 pretreatment (Fig.7a). This also identified prac- 
tolol and oxprenolol, both of which led to a suppressed DMR by 
U0126 (Supplementary Fig.Sll). Second, we examined the impacts 
of a Casein kinase (CK2) inhibitor TBB on the drug-induced DMR 
(Fig. 7b), wherein TBB was used to precondition the cells with sup- 
pressed CK2 activity. Distinct drugs exhibited different sensitivities 
to the TBB treatment (Supplementary Fig.S12). Third, we compared 
a drug-induced DMR with its whole cell cAMP signal (Fig.8). This 
identified DMR-biased drugs including practolol, alprenolol, pindo- 
lol, labetalol, acebutolol, and dopamine. Dopamine was previously 
shown to be moderately potent to trigger DMR signal in A431 cells 29 . 

Discussion 

In contrast to conventional molecular assays that are biased towards 
a single pathway and/or a single molecule, the integrative readout of 
DMR assays allows assessing drug pharmacology with wide pathway 
coverage. The imminent high temporal resolution makes DMR 
assays possible to quantify drug pharmacology at different time 
domains. The non-invasive measure enables DMR assays to probe 
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Figure 3 | Heat maps of P-drugs. (a) Clinical indications of all drugs were used as a basis for similarity analysis, (b) Chemical structures of all drugs were 
used as a basis for similarity analysis. 
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Figure 4 | Correlation analysis identifies the origin of similarity among 
distinct f-drugs based on potency, mechanisms of activation and 
deactivation, (a) Scatter plots between the |3-drug-induced DMR in native 
cells and the EPI-induced DMR after pretreated with drugs for lhr (black 
dots) or 5hr (pink dots); and (b) scatter plots between the (3-drug-induced 
DMR in cells after pretreatment with EPI for lhr and 5hr. 
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Figure 5 | Correlation analysis identifies the origin of similarity among 
distinct 0-drugs based on potency, mechanisms of activation and 
deactivation, (a) Scatter plots between the |3-drug-induced DMR without 
and with the pretreatment with PRO for lhr; (b) scatter plots between the 
[3-drug-induced DMR in native cells and the PRO-induced DMR in the 
cells after 5hr pretreatment with (3-drugs. 



drug pharmacology under a wide range of conditions. The DMR 
assays are robust and reproducible (Supplementary Fig.S13). 
Combining DMR assays with similarity analysis can differentiate 
drug pharmacology with relatively high resolutions. We implanted 
the iPOT to investigate the pharmacology of clinically available p- AR 
drugs. Using twelve DMR assays in a cell line endogenously expres- 
sing the p 2 - AR under distinct chemical environments, we found that 
clinically approved P-AR drugs are divergent in their potency, 
mechanisms of activation and inactivation, and pathway biased 
activity. Although only a first step, our results are encouraging. 
The most notable finding is that the iPOT profiling of P-AR drugs 
established an apparently effective link between in vitro results and 
in vivo indications. iPOT profiles of drugs in disease relevant cells, 
together with global analysis of drugs against the family of adrenergic 
receptors, would further strengthen the correlation between in vitro 
and in vivo results. 

Since DMR is an integrated measure of receptor signaling, specific 
biological events that underlie the similarity and differences among 
distinct drugs are still largely unknown 30 . The linkage between a 
specific DMR parameter and a specific clinical feature is also largely 
unknown at the present time. Further delineation of receptor biology 
and drug pharmacology, together with optimization of algorithm for 
similarity analysis, would be necessary to validate the link between 
the iPOT testing results and in vivo indications, and to make it 
possible to rank and select drugs within a class for in vivo testing 
during lead selection process. Nonetheless, the high content assess- 



ment of drug pharmacology with label-free iPOT assays elucidates 
the complex and multifaceted efficacy of GPCR ligands, thus offering 
a practical process to reposition drugs for new indications 31,32 and to 
prioritize lead molecules. 

Methods 

Cell lines and reagents. Human epidermoid carcinoma A431 was obtained from 
American Type Cell Culture. The cells were grown in Dulbecco's modified Eagle's 
medium (DMEM) supplemented with 10% fetal bovine serum (FBS), 4.5g/liter 
glucose, 2 mM glutamine, and antibiotics. The native cells were passed when 
approaching —90% confluence with trypin/EDTA (etthylenediaminetetraacetic acid) 
to provide new maintenance cultures on T-75 flasks and experimental cultures on the 
biosensor microplates. All drug molecules were purchased from commercial sources 
(Supplementary Table S2). All drugs were dissolved in DMSO and were diluted 
directly into the assay buffer (IX Hanks' balanced salt buffer, 20mM Hepes, pH 7.1; 
HBSS) to the indicated concentrations. Epic® 384-well tissue culture treated 
biosensor microplates were obtained from Corning Incorporated (Corning, NY, 
USA). 

DMR assays. Epic® system (Corning Inc.), a wavelength interrogation reader system 
tailored for resonant waveguide grating (RWG) biosensors in microtiter plates, was 
used for all DMR assays. This system consists of a temperature-control unit (26°C), 
an optical detection unit, and an on-board liquid handling unit with robotics. The 
detection unit is centered on integrated fiber optics, and enables kinetic measures of 
cellular responses with a time interval of ~15sec. 20,000 cells per well were directly 
seeded and cultured overnight under serum rich medium, followed by starvation 
overnight using serum depleted medium. After washed twice, the cells were 
maintained with HBSS and further incubated inside the system for lhr. A 2-min 
baseline was then established. Immediate after the compound additions using the 
onboard liquid handler, the cellular responses were recorded. All studies were carried 
out with at least three replicates. 
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Figure 6 | Correlation analysis identifies the origin of similarity among 
distinct f-drugs based on pathway biased agonism. (a) Scatter plot 
between the p-drug-induced DMR in native cells and the FSK DMR in the 
presence of distinct (3-drugs; (b) Scatter plot between the P-drug-induced 
DMR in native cells and in the cells after pretreatment with CTx. 



Drug screening. We screened all 36 clinically available (3-AR drugs and drug isomers 
against 12 distinct assay conditions using kinetic DMR assays. Drug concentrations 
were chosen to achieve maximum efficacy and to be amenable for high throughput 
screening (Supplementary Table S3). A431 cells were preconditioned with probe 
molecules to achieve a wide range of chemical environments, which, in turn, manifest 
the specificity, relative potency and efficacy, and modes of action of the drugs. The 
chemical environments also enable mapping the cell systems -based, ligand- directed 
and kinetics-dependent biased agonism of the drugs. Specifically, cells were 
pretreated with DMSO (the positive control), 5nM epinephrine, luM propranolol, 
400ng/ml cholera toxin, lOOng/ml pertussis toxin, U0126, 4,5,6,7- 
tetrabromobenzotriazole, forskolin, and the drug molecules for a specific period of 
time (Table 1). Afterwards, cells were stimulated with a (3 2 -hgand (epinephrine, 
propranolol, or the drug), whose responses were recorded in real time and used for 
similarity and correlation analysis. 

Whole cAMP assays. Cells were plated in 384well tissue culture treated plates (BD 
Bioscience, Cat# 354660) with a seeding density of 20000 cells per well. Cells were 
cultured in the serum rich medium overnight. The next day cAMP-Glo assay was 
performed according to manufacturer's instruction (Promega, Cat#V1502). Cells 
were incubated with 10 uM compounds in induction buffer for 30 minutes before 
adding lysis buffer. Luminescence was measured using Tecan Safirell reader. 

Data visualization and clustering. For each DMR the responses at the four distinct 
time points (Fig.lc) were extracted for reduction of the time dimensions. All time 
points refer to the stimulation duration after renormalized the responses starting 
from t 0 (the time when the compound was added) (Fig.lc). For visualization purpose 
the responses were color coded to illustrate relative differences in DMR signal 
strength. The total 48 dimensions, 4 for each condition, were thus obtained to rewrite 
the DMR pharmacology of a drug. In the drug matrix (Fig. 2) each column represents 
one DMR response at a particular time in a specific assay condition, and each row 
represents one drug with the original approved indications, as documented in 
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Figure 7 | Correlation analysis identifies the origin of similarity among 
distinct f-drugs based on pathway biased agonism. Scatter plots between 
the p-drug-induced DMR in native cells and in the cells after pretreatment 
with PTx (black dots) or U0126 (pink dots) (a), and TBB (b). 



DrugBank (Http:// www.drugbank.ca). Every row and column carries equal weight. 
The Ward hierarchical clustering algorithm and euclidean distance metrics (http:// 
www.eisenlab.org/eisen/) were used for clustering the results. DMSO in the vehicle, a 
concentration that equals to those for all drugs, is also included as a negative control. 
For replicates, we used statistical differences using the four chosen time points to 
calculate the variability among replicates. All replicates have to pass the 2 X coefficient 
of variation test before being included in analysis. The drugs whose DMR responses 
failed the test are re- screened. For clinical indications, each drug was written using a 
binary numeral system, 1 for a clinical indication as documented in DrugBank and 0 
for none. For chemical structures, all drugs were analyzed using ChemMine online 
software, which employs atom pairs as structural descriptors and Tanimoto 
coefficient as a similarity measure. Hierarchical clustering was used for clustering 
based on both clinical features and chemical structures. 
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